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Abstract — Heat  transfer  is  one  of  a  number  of 
indispensable  tools  in  studying  of  ICEs,  due  to  its 
influence  on  decisive  parameters  of  operation  like 
temperature  and  pressure  inside  the  cylinder.  It  is  safe  to 
say  that  modeling  of  the  engine  heat  transfer  is  among  the 
most  complex  problems  for  engineers.  Application  of 
numerical  methods  to  predict  the  heat  transfer  in  a 
cylinder  of  reciprocating  ICEs  is  a  process  of  high 
importance,  which  was  recognized  from  the  earliest 
stages  of  their  development.  This  is  done  to  examine 
performance  optimization  and  design  improvement  in 
order  to  meet  nowadays  demands  exhibited  on  the 
engines.  This  present  study  focuses  on  a  3-D  transient 
state  temperature  distribution  analysis  on  a  gasoline 
engine  model  via  formulated  of  models,  simulating  using 
EEM  in-built  in  the  COMSOL  Multiphysics  software  4.3a 
to  determine  the  temperature  distribution  and  gradient  of 
the  engine  cylinder  head  model.  The  number  of  degrees  of 
freedom  solved  for  were  32685  in  383  s  (12  minutes,  20 
seconds)  in  the  mesh  optimization.  From  the  result 
obtained  it  was  discovered  that  the  heat  transfer  in  the 
combustion  chamber  of  the  ICE  varies  with  time.  Thus,  it 
took  the  engine  10  minutes  to  complete  a  cycle  vis-a-vis 
transfer  of  heat  after  combustion  and  that  the  heat 
transfer  starts  after  30  seconds  of  combustion.  In 
addition,  the  temperature  of  the  cylinder  dropped  from 
1273.2  K  to  301  K  over  a  period  of  10  minute. 

Keywords —  Comsol,  Heat  transfer,  ICE ,  Simulation , 
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I.  INTRODUCTION 

One  of  the  greatest  contributors  to  engine  energy  losses  is 
heat  loss.  Almost  one  third  of  the  heat  released  during 
combustion  is  lost  to  cooling  media  and  exhaust  gases, 
thus  dominating  the  Internal  Combustion  Engine  (ICE) 
energy  balance  (Dahlstrom,  2016).  The  heat  transfer  from 
combustion  gases  to  the  coolant  in  reciprocating  ICEs 
varies  between  25  -  35%  of  the  total  energy  released  by 
the  mixture  of  fuel  and  combustion  air.  Nearly  half  of  the 


heat  passes  through  the  cylinder  walls  and  most  of  the 
remaining  heat  is  transferred  to  the  coolant  in  the  cylinder 
head  with  the  most  intensive  heat  transfer  near  the 
exhaust  valve  seats  (Oleg,  2013). 

During  engine  operation,  there  are  energy  losses  via  heat 
transfer  to  the  cylinder  walls  during  the  latter  stages  of  the 
compression  stroke  and  the  combustion-expansion  stroke. 
During  the  compression  process,  the  charge  temperature 
rises  above  the  wall  temperature  and  thus  heat  transfer 
occurs  from  the  in-cylinder  gases  to  the  cylinder  walls.  In 
the  combustion-expansion  stroke,  the  heat  transfer  rate  to 
the  cylinder  walls  is  the  highest  (Wei  et  al .,  2016). 


intake 


compression 


power 


exhaust 


Fig.  1:  A  four-stroke  ignition  cycle  ( Antonio ,  2008) 


Modeling  of  the  heat  transfer  is  usually  considered  as  a 
part  of  whole  engine  simulation  including  combustion 
modeling,  emission  and  soot  formation  modeling,  fuel 
spray  modeling  etc.,  and  can  serve  as  a  prerequisite  of 
performance  optimization  and  design  improvement  in 
order  to  meet  nowadays  demands  added  on  the  engines. 
There  are  two  basic  types  of  models  that  are  in  use  for 
modeling  of  the  heat  transfer:  thermodynamic  and  multi¬ 
dimensional  models.  The  simplest  type  of  thermodynamic 
models  are  zero  dimensional  models  or  global  models, 
which  suits  for  estimating  of  the  general  parameters  by 
use  of  semi  empiric  correlations.  Another  kind  of 
thermodynamic  models  is  quasi  dimensional  models. 
Currently,  researchers  tend  to  move  towards  more 
comprehensive  and  accurate  models  describing  the 
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performance  of  the  engine,  although,  in  some  cases 
proved  to  be  impractical  for  their  complexity.  Multi¬ 
dimensional  models  are  increasingly  gaining  more 
attention  due  to  their  ability  to  compute  the  gas  motion  by 
numerical  solving  of  the  differential  equations 
representing  the  laws  of  mass,  energy  and  momentum 
conservation,  thus,  providing  a  greater  deal  of  precision 
(Oleg,  2013). 

Global  approach  was  employed  by  Ziarati  (1990)  in  1996 
to  study  the  delay  period  and  droplet  penetration.  Thanks 
to  the  developed  model  which  used  the  Annand’s 
correlation  to  describe  the  heat  transfer  in  the 
combustion  chamber  and  was  tested  against  the 
experimental  results  of  a  direct  injection  diesel  engine 
(Ricardo,  Atlas)  for  a  number  of  nozzle  configuration, 
plunger  size  and  engine  speeds  and  loads.  Zarati  managed 
to  obtain  good  agreement  of  cylinder  pressure  and 
temperature  with  experimental  results 
However,  many  attempts  were  made  to  study  and 
understand  the  causes  of  these  heat  losses  and  what 
amendments  could  be  done,  mainly  through  simulations 
and  experimental  work  (Dahl  strom,  2016).  According  to 
Vijayakumaran  et  al.  (2014),  any  physical  fluid  flow 
problem  can  be  solved  either  experimentally  or 
numerically.  In  the  recent  years  there  were  tremendous 
improvement  in  the  field  of  numerical  techniques  which 
made  a  great  impact  on  the  evaluation  of  complex  flow 
problems  and  achieving  their  solution. 

Xu  and  Cho  (2016)  reported  that  In-cylinder  flow 
simulation  experiment  works  have  been  carried  out  using 
a  hot  wire  velocimetry  or  laser  Doppler  velocimetry 
(LDV)  to  measure  the  velocity  field.  This  test  is  very 
difficult  because  the  characteristics  measured  in  the 
reciprocating  engine  cylinder  flow  are  a  highly  complex 
three-dimensional  turbulence  and  instability.  The  paper 
therefore  concluded  that  the  numerical  approach  could  be 
changed  by  developing  a  CFD  model  for  the  in-cylinder 
flow  predictions. 

Kumiawan  and  Abdullah  (2008)  performed  the  CFD 
simulation  for  the  combustion  process  in  their  analysis 
using  the  moving  mesh  and  boundary  algorithm  where 
every  event  represents  the  different  mesh  and  boundary 
geometries  for  every  different  crank  angle  in  each  step  of 
engine  cycle.  Thus,  in  order  to  perform  CFD  simulation 
for  internal  combustion  process,  the  analysis  was  carried 
out  by  using  transient,  moving  meshes  and  boundaries, 
high  compressible  Reynolds  number,  high  turbulence 
intensity,  momentum,  heat  and  mass  transfer  with 
complex  geometries  model  and  chemical -thermal 
dependent  as  well. 

As  reported  by  Meman  et  al.  (2016),  Cheng  x.  et  al. 
(2015)  analyzed  on  the  cylinder  head  of  the  heavy  diesel 
engine  by  performing  3 -dimensional  CFD  analysis  with 
the  cylinder  combustion  and  an  engine  cooling  jacket, 
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stress  and  strain  analysis  under  thermal -mechanical 
coupling  conditions  and  high  cycle  fatigue  analysis  and 
also  they  had  measured  temperature  distribution  on  the 
engine  head  and  predicted  the  danger  region  through 
calculation  matches  the  crack  in  actual  experiment.  New 
thermo-mechanical  analysis  of  cylinder  head  was 
investigated  by  Fadaei  et  al.  (2011)  in  this  paper  the  3D 
model  was  created  by  using  Solid  Works  2008.  The  mesh 
and  analysis  of  cylinder  head  were  constructed  in 
ANSYS. 

Meman  et  al.  (2016)  also  reported  the  work  of  Baox  in 
Zhao  et  al.  (2012)  and  the  paper  evaluated  the  problem  of 
super  heat  generation  in  cylinder  head.  The  simulation  of 
heat  transfer  was  taken  out  for  the  combustion  chamber 
by  help  of  the  coupled  finite  element  method  (FEM). 
Hence,  the  temperature  distribution  at  cylinder  blocks, 
cylinder,  cylinder  head,  cylinder  gasket,  valve  seat  and 
valve  guide  were  measured. 

Oke  et  al.  (2015)  confronted  the  problem  of  heat  transfer 
in  spark  ignition  engine  by  developing  a  model  for  2- 
stroke  spark  ignition  engine  and  used  Comsol 
Multiphysics  to  predict  the  temperature  distributions  and 
heat  transfer  on  the  combustion  chamber. 

Despite  all  these  studies  regarding  the  effect  of  heat 
transfer  on  the  combustion  chamber  walls  of  Internal 
Combustion  Engines  (ICE),  there  are  still  areas  which 
have  yet  to  be  examined.  This  work  concerns  heat  transfer 
principles  in  four  strokes  ICE  from  modeling  and 
simulation  standpoint. 

II.  METHODOLOGY 

2.1  Modeling  and  Simulation 

2.1.1  Model  Formulation  -  One  Dimensional 
Solution  to  the  Problem 

Domain  Equation 

Starting  with  the  time  dependent  heat  conduction 
equation 

pCpJ  +  v.(-/cvr)  =  (?  (i) 

in  Cartesian  form 

pcp  dt  +  dA  kx  dx )  ~  Q 

To  avoid  division  by  zero  at  x=0,  which  creates  numerical 
problems,  multiply  equation  (2)  by  x2 


Using  dimensionless  coordinate  v  by  scaling  the  equation 
provides  the  option  to  quickly  change  the  wall  thickness 
without  changing  or  parameterizing  the  geometry  size. 
Introducing  the  dimensionless  coordinate, 
x  d  Id 
wx  ’  dx  wx  dx 
and  substituting  in  equation  (iii)  leads  to 

.9  7  dT  1  9  (  .9  9  .9  9 

*  W  =  1  W*Q 
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-2  dT 
x  Pcv  iF 
On  the  domain 

x=  0 


,  d  (  kxx2  d T\  a2/->  x 

+  ^ - ^  )=XQ  (!V) 

dx  \  wx  dxj 


1=1 


I - 1* 

Comsol  implementation 

To  implement  equation  (iv)  and  the  boundary  conditions 
of  this  problem,  an  ID  time -dependent  equation  version 
of  the  General  Form  PDE  interface  given  in  equation  (v) 
was  used: 


d  u  .  , 

ea  dt2  da 


du 

dt 


+  v.r 


— nr  —  G  + 

0  =  wY 


fdWxY 

(id  p 


inQ 

on  df2 
on  on 


(v) 


Comparing  equation  (iv)  to 

d2u  du 

eaatr  +  daat  +v-r  =  F 
ea  =  0 
da  =  x2pcp 


r  =  -■ 


kxx2  dT 
wr  dx 


Q  =  0 

Boundary  and  initial  conditions 

Because  of  symmetry  about  x=0,  there  is  no  flux  through 
this  point  meaning 

dT 

d2  =  ° 

At  the  surface  x=l,  a  convective  heating  and  cooling 
expression  with  a  heat  transfer  coefficient  hs  [W/(m2.K)], 
for  the  influx  of  heat  (W/m2). 

Qin  =  K  (Text  ~  *0 

Since  the  flux, 


r  =  - 


kxx2 dT 
wr  dx 


and 


Therefore, 


dT 

~kx  —  hs  (Text  ~  T') 


9  —  hs(Text  T) 

2.1.2  COMSOL  Multiphysics  4.3a 

COMSOL  Multiphysics  4.3a  software  package  which  is  a 
powerful  finite  element  method  (FEM)  and  partial 
differential  equation  (PDE)  solution  engine  was  used  for 
the  analysis.  It  has  ten  (10)  add-on  modules  that  expand 
the  capabilities  in  the  following  application  areas  (Alabi 
et  al. ,  2015):  AC/DC,  Acoustics,  Chemical  Species  and 
Transport,  Electrochemistry,  Fluid  flow,  Heat  transfer, 
Plasma,  Radio  Frequency  (RF),  Structural  Mechanics, 
and  Mathematics  Modules. 


Fig.  2  below  is  the  schematic  representation  of  Comsol 
Multiphysics  processes  or  steps  as  applicable  to  this 
work. 


Fig.  2:  Comsol  Schematic  Representation  of  the  Study 


1.  Model  Pre-processing 

The  Structural  Mechanics  add-on  module  which  has 
a  Thermal-Strength  in-built  model  was  used  for  the 
simulation  and  analysis. 

pCp5  =  v.(/cvr)  (vi) 


Fig. 3:  Geometry  model  window  in  Inventor 
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Fig.  4:  Mesh  representation  of  the  study 

2.  Model-Processing 

Mesh  Generation:  The  mesh  chosen  was  Free 
Tetrahedral.  The  Completed  mesh  consists  of  larger 
number  of  elements  and  degrees  of  freedoms  in  the 
best  possible  period. 

3.  Model-Post-processing  and  Visualization 

In  post  processing  mode  different  plots  types  were 
obtained  from  sets  of  selected  parameters.  The 
solution  variables  were  later  visualized  using  the  post 
processing  utilities  for  their  derivative,  and  the  space 
coordinates.  Many  frequently  used  expressions  were 
predefined  as  application  mode  variables,  directly 
available  from  lists  in  the  Plot  Parameters. 

III.  RESULTS  AND  DISCUSSION 

The  results  are  presented  as  shown  in  the  figures  below. 


Fig.  5:  Full-cylinder  head  Model  at  initial  combustion 
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Fig.  6:  Model  after  30  seconds  of  combustion 
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Fig.  5  shows  the  full  cylinder  model.  The  cylinder 
reached  its  maximum  temperature  of  1273.2K  at  the 
earliest  stage  of  combustion  which  is  when  the 
combustion  just  occurs.  After  0.5  minutes  heat  transfer 
process  began  and  traces  of  red  hot  start  disappearing  and 
this  could  be  seen  in  fig.  6. 

From  fig.  7  more  disappearance  of  heat  was  pronounced 
and  After  9  minutes  the  traces  of  heat  as  shown  in  fig.  8 
has  completely  disappear.  Fig.  9  shows  the  time 
dependent  solver  while  fig.  10  is  the  result  obtained  from 
the  1-D  heat  transfer  model,  which  is  the  plot  of 
temperature  against  time.  It  was  observed  that  the 
temperature  of  the  cylinder  dropped  from  1273.2  K  to  301 
K  over  a  period  of  10  minute.  This  is  an  indication  that  a 
complete  cycle  took  10  minute  to  transfer  heat  after 
combustion,  and  that  the  heat  transfer  starts  after  30 
seconds  of  combustion. 


Fig.  7:  Heat  Transfer  after  4  minutes 
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Fig.  9:  Time  Dependent  Solver 


Fig.  10:  Plot  of  temperature  against  time 
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IV.  CONCLUSION 

The  energy  losses  resulting  from  the  transfer  of  heat  to 
the  cylinder  walls  is  mainly  through  the  piston  crown, 
cylinder  liner,  combustion  chamber,  and  valves.  A 
decrease  in  the  surface  temperature  of  these  components 
leads  to  a  reduction  of  the  transfer  of  heat  to  the  engine 
coolant  and  thus  reduces  energy  losses  (Wei  et  al. ,  2016). 
Thus,  the  heat  transfer  in  the  combustion  chamber  of  the 
ICE  varies  with  respect  to  time.  The  obtained  results 
showed  that  it  took  the  engine  10  minutes  to  complete  a 
cycle  vis-a-vis  transfer  of  heat  after  combustion,  and  that 
the  heat  transfer  starts  after  30  seconds  of  combustion.  In 
addition,  the  temperature  of  the  cylinder  dropped  from 
1273.2  K  to  301  K  over  a  period  of  10  minute. 

An  advancement  in  technology  leads  to  this  research,  for 
improvement  and  longevity  in  four  stroke  engine 
performance.  The  results  obtained  can  be  validated  using 
other  available  software  packages  like  ANSYS, 
ABAQUS,  NASTRAN,  MATLAB,  J-ANALYZER  used 
in  simulating  the  heat  transfer  on  IC  Engines. 
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